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(^— N ■ We describe a quantum electromechanical system(QEMS) comprising a single quantum dot har- 

*vj ' monically bound between two electrodes and facilitating a tunneling current between them. An 

example of such a system is a fuUerene molecule between two metal electrodes [Park et al. Nature, 

C . 407, 57 (2000)]. The description is based on a quantum master equation for the density operator 

, ^ ' of the electronic and vibrational degrees of freedom and thus incorporates the dynamics of both 

diagonal (population) and off diagonal (coherence) terms. We derive coupled equations of motion 

for the electron occupation number of the dot and the vibrational degrees of freedom, including 

damping of the vibration and thermo-mechanical noise. This dynamical description is related to 

observable features of the system including the stationary current as a function of bias voltage. 
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PACS numbers: 72.70.+m,73.23.-b,73.63.Kv,62.25.+g,61.46. 



I. INTRODUCTION 



A quantum electromechanical system (QEMS) is a submicron electromechanical device fabricated through state- 
of-the-art nanofabrication^ . Typically, such devices comprise a mechanical oscillator (a singly or doubly clamped 
cantilever) with surface wires patterned through shadow mask metal evaporation. The wires can be used to drive the 
i-rt ! mechanical system by carrying an AC current in an external static magnetic field. Surface wires can also be used as 
P5 > motion transducers through induced EMFs as the substrate oscillates in the external magnetic field. Alternatively 
O the mechanical resonators can form an active part of a single electron transducer, such as one plate of a capacitively 
^, t coupled single electron transistoriM,. These devices have been proposed as sensitive weak force probes with the potential 
to achieve single spin detectionai^. However they are of considerable interest in their own right as nano-fabricated 
— I , mechanical resonators capable of exhibiting quantum noise features, such as squeezing and entanglement* . 
^ ' In order to observe quantum noise in a QEMS device we must recognize that these devices are open quantum 

vIN [ systems and explicitly describe the interactions between the device and a number of thermal reservoirs. This is the 
primary objective of this paper. There are several factors that determine whether a system operates in the quantum 
or classical regime. When the system consists only of an oscillator coupled to a bath the oscillator quantum of 
^^ ,. energy should be greater than the thermo-mechanical excitation of the system; hujo > ksT where Wo is the resonant 
^+ ' frequency of the QEMS oscillator and T is the temperature of the thermal mechanical bath in equilibrium with the 
(T^ , oscillator. At a temperature of 10 milliKelvin, this implies an oscillator frequency of the order of GHz or greater, 
^i Recently Huang et al. reported the operation of a GHz mechanical oscillator—. A very different approach to achieving 
2 . cL high mechanical frequency was the fuUerene molecular system of Park et al.S,, and it is this system which we take 
as the prototype for our theoretical description. Previous work on the micro-mechanical degrees of freedom coupled 
' ' to mesoscopic conductorsSiifliiiiiS, indicate that transport of carriers between source and drain can act as a damping 
~^ I reservoir, even in the absence of an other explicit mechanisnr for mechanical damping into a thermal reservoir. This 
Q ' is also predicted by the theory we present for a particular bias condition. As is well known, dissipation can restore 
^ , semiclassical behavior. Transport induced damping can also achieve this result. 

The model we describe, Fig. ^ consists of a single quantum dot coupled via tunnel junctions to two reservoirs, 

the source and the drain. We will assume that the coulomb blockade permits only one quasi-bound single electron 

^ • state on the dot which participates in the tunneling between the source and the drain. We will ignore spin, as the 

J-j [ source and drain are not spin polarized, and there is no external magnetic field. A gate voltage controls the energy 

. . . 1 of this quasi-bound state with respect to the Fermi energy in the source. The quantum dot can oscillate around 

an equilibrium position mid way between the source and the drain contacts due to weak restoring forces. When 

an electron tunnels onto the dot an electrostatic force is exerted on the dot shifting its equilibrium position. In 

essence this is a quantum dot single electron transistor. In the experiment of Park et al.^, the quantum dot was a 

single fuUerene molecule weakly bound by van der Walls interactions between the molecule and the electrodes. The 

dependence of the conductance on gate voltage was found to exhibit features attributed to transitions between the 

quantized vibrational levels of the mechanical oscillations of the molecule. 

Boese and Schoeller^'^ have recently given a theoretical description of the conductance features of this system. A 
more detailed analysis using similar techniques was given by Aji et al^. Our objective is to extend these models to 
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FIG. 1: Schematic representation of tunneling between a source and a drain through a quantum dot. The dot is harmonically 
bound and vibrational motion can be excited as electrons tunnel through the system. 

provide a full master equation description of the irreversible dynamics, including quantum correlation between the 
mechanical and electronic degrees of freedom. We wish to go beyond a rate equation description so as to be able 
to include coherent quantum effects which arise, for example, when the mechanical degree of freedom is subject to 
coherent driving. 

II. THE MODEL 

We will assume that the center of mass of the dot is bound in a harmonic potential with resonant frequency ojo- 
This vibrational degree of freedom is then described by a displacement operator x which can be written in terms of 
annihilation and creation operators a, a^ as 
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The electronic single quasi-bound state on the dot is described by Fermi annihilation and creation operators c, c' , 
which satisfy the anti commutation relation cc^ + c^c — 1. 
The Hamiltonian of the system can then be written as, 
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where n = c^c is the excess electron number operator on the dot. 

The first term of the Hamiltonian describes the free energy for the island. A particular gate voltage Vg, with a 
corresponding huji = IbmeV, for the island is chosen for calculation. Uc is the Coulomb charge energy which is the 
energy that is required to add an electron when there is already one electron occupying the island. We will assume 
this energy is large enough so that no more than one electron occupies the island at any time. This is a Coulomb 
Blockade effect. The charging energy of the fuUerene molecule transistor has been observed by Park et al. to be 
larger than 270 meV which is two orders of magnitude larger than the vibrational quantum of energy hujo- The free 
Hamiltonian for the oscillator is described in ternr ^. The Park et al. experiment gives the value huJo = 5 meV, 
corresponding to a THz oscillator. The electrostatic energy of electrons in the source and drain reservoirs is written 
as term Q. Term (O is the coupling between the oscillator and charge while term (jSJ represents the source-island 
tunnel coupling and the drain- island tunnel coupling. The last term, ((TJ, describes the coupling between the oscillator 
and the thcrmo-mechanical bath responsible for damping and thermal noise in the mechanical system in the rotating 
wave approximation^^. This is an additional source of damping to that which can arise due to the transport process 



itself (see below). We include it in order to bound the motion under certain bias conditions. A possible physical 
origin of this source of dissipation will be discussed after the derivation of the master equation. 

We have neglected the position dependence of the tunneling rate onto and off the island. This is equivalent to 
assuming that the distance, d, between the electrodes and the equilibrium position of the uncharged quantum dot, is 
much larger than the rms position fluctuations in the ground state of the oscillator. There are important situations 
where this approximation cannot be made, for example in the so called 'charge shuttle' systemsii. 

A primary difficulty in analyzing the quantum dynamics of this open system is the presence of different time scales 
associated with the oscillator, the tunneling events and the coupling between the oscillator and electronic degrees of 
degree due to the electrostatic potential, term (4). The standard approach would be to move to an interaction picture 
for the oscillator and the electronic degrees of freedom. However this would make the electrostatic coupling energy 
time dependent, and rapidly oscillating. Were we to approximate this with the secular terms stemming from a Dyson 
expansion of the Hamiltonian, the resulting effective coupling between the oscillator and the electron occupation of 
the dot simply shifts the free energy of the dot and no excitation of the mechanical motion can occur. 

To avoid this problem we eliminate the coupling term of the oscillator and charge by doing a canonical transformation 
with unitary representation U — e'^ where, 

s = -A(a'l' - a)h (8) 



with 
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This unitary transformation gives a conditional displacement of the oscillator conditional on the electronic occupation 
the dot. One might call this a displacement picture. 

This derivation follows the approach of MaharJ>£. The motivation behind this is as follows. The electrostatic 
interaction, term (4), displaces the equilibrium position of the oscillator so that the average value of the oscillator 
amplitude in the ground state becomes 

(a) = A, (10) 

We can shift this back to the origin by a phase-space displacement 

a = e^ae"^ = a + An. (11) 

This unitary transformation gives a conditional displacement of the oscillator, conditional on the electronic occupation 
of the dot. Applying U to the Fermi operator c gives 

c = ce^(''*-''). (12) 

The Schrodinger equation for the displaced state, p — e'^pe^^, then takes the form 



where the transformed Hamiltonian is 
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We will now work exclusively in this displacement picture. 

To derive a master equation for the dot, we first transform to an interaction picture in the usual way to give the 
Hamiltonian 

Hi = ^rsfe(afec^e*('^^~''~'^^^)*e~^('''""'°'^'^'=""'°') 
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where i] = y^ /{huo) = X-^- 

At this point we might wish to trace out the phonon bath, however we will postpone this for a closer look at the 
tunneling Hamiltonian at the individual phonon level. We use the exponential approximation e^ = l + a; + a;^/(2!) + - • • 
, when X is small for the term e^'^'^ '^ ° ~°^ ° ^ . We expect an expansion to second order in A to give an adequate 
description of transport, in that at least one step in the current vs. bias voltage curve is seen due to phonon mediated 
tunneling. In the experiment of Park et al., A was less than unity, but not very small. Strong coupling between the 
electronic and vibrational degrees of freedom (large A) will give multi-phonon tunneling events, and corresponding 
multiple steps in the current vs. bias voltage curves. The Hamiltonian can then be written as 
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The terms of zero order in A describe bare tunneling through the system and do not cause excitations of the vibrational 
degree of freedom. The terms linear in A correspond to the exchange of one vibrational quantum, or phonon. The 
terms quadratic in A correspond to tunneling with the exchange of two vibrational quanta. Higher order terms could 
obviously be included at considerable computational expense. We will proceed to derive the master equation up to 
quadratic order in A. 

III. MASTER EQUATION 

Our objective here is to find an evolution equation of the joint density operator for the electronic and vibrational 
degrees of freedom. We will use standard methods based on the Born and Markov approximationpi^. In order to 
indicate where these approximations occur, we will sketch some of the key elements of the derivation in what follows. 
The Born approximation assumes that the coupling between the leads and the local system is weak and thus second 
order perturbation theory will suffice to describe this interaction; 

P ^ ^ I di'Tr[Hi(t),[Hi(t'),R]], (17) 

n Jo 

where R is the joint density matrix for the vibrational and electronic degrees of freedom of the local system and the 
reservoirs. 

At this point we would like to trace out the electronic degrees of freedom for the source and drain. We will assume 
that the states of the source and drain remain in local thermodynamic equilibrium at temperature T. This is part of 
the Markov approximation. Its validity requires that any correlation that develops between the electrons in the leads 
and the local system, as a result of the tunneling interaction, is rapidly damped to zero on time scales relevant for 
the system dynamics. We need the following moments: 

Tr[a|.afcp] = fsk, Tr[6^6fep] = fok, 

Tr[afea^p] = 1- fsk, T:r[bkb\p\ = 1 - fok- 



where fsk = f{Esk) is the Fermi function describing the average occupation number in the source and similarly 
fok = f{EDk), for the drain. The Fermi function has an implicit dependence on the temperature, T, of the electronic 
system. 

The next step is to convert the sum over modes to a frequency-space integral: 

y^fsk\Tsk?^ dujg{uj)fD{uj)\Ts{io)\^, (18) 

k -^0 

where iTskl"^ — Tgf.Tsk and g{u!) is the density of states. Evaluating the time integral, we use, 

/•OO 

dre^''^ = 7r<5(e) ± iPV{l/e), (19) 

where t = t — t' and the imaginary term is ignored. 

Using these methods, we can combine the terms for the source and drain as the left and right tunneling rates, 7l 
and 7/f respectively 



doj g{uj)\Ts{oj)Y5(u:Q) = 7l(wo). (20) 





In the same way, we can define 



7-Ll = iLiflLOi -T]- flL), 

fiL = fi'hoJi - ?7 - fiL), 

fsL = f{huji -r] + hujo- ^iL)■ 

and similarly for "fRi,"fR2"fR3, fiR, f2R, fsR replacing fiL with iir and / being the Fermi functions which have a 
dependence on the bias voltage (through the chemical potential) and also on the phonon energy hujo. As the bias 
voltage is increased from zero, the first Fermi function to be significantly different from zero is /2L followed by /i^ 
and then f^^ ■ This stepwise behavior will be important in understanding the dependence of the stationary current 
as a function of bias voltage. 

The master equation in the canonical transformed picture to the second order in A may be written as 
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where the notation T> is defined for arbitrary operators X and Y as 

V[X]Y = J[X]Y-A[Y] 



= XYX^ - Ux^XY + YX^X), (22) 



and 



^P('^°) = ^huJksT _ I ■ (23) 

We have included in this model an explicit damping process of the oscillators motion at rate k into a thermal oscillator 
bath with mean excitation rip and T is the effective temperature of reservoir responsible for this damping process. 
A possible physical origin for this kind of damping could be as follows. Thermal fluctuations in the metal contacts 
of the source and drain cause fluctuations in position of the center of the trapping potential confining the molecule, 
that is to say small, fluctuating linear forces act on the molecule. For a harmonic trap, this appears to the oscillator 
as a thermal bath. However we expect such a mechanism to be very weak. This fact, together with the very large 
frequency of the oscillator, justifles our use of the quantum optical master equation (as opposed to the Brownian 
motion master equation) to describe this source of dissipation-^^. The thermo-mechanical and electronic temperatures 
are not necessarily the same, although we will generally assume this to be the case. 

Setting A = we recover the standard master equation for a single quantum dot coupled to leadsi^. The superop- 
erator ^[c^] adds one electron to the dot. Terms containing this super operator describe a conditional Poisson event 
in which an electron tunnels onto the dot. The electron can enter from the source, with probability per unit time 
of 7li/il(ccT), or it can enter from the drain, with probability per unit time 7_ri/i_r (cc^). Likewise the term 2?[c] 
describes an electron leaving the dot, again via tunneling into the source (terms proportional to ^li) or the drain 
(terms proportional to 7_ri). When A 7^ there are additional terms describing phonon mediated tunneling events 
onto and off the dot. Any term proportional to 7^^, i = 1,2,3. describes an electron tunneling out of, or into, the 
source, while any term proportional to 7^^, i = 1, 2, 3 describes an electron tunneling out of, or into, the drain. 

The average currents through the left junction (source lead-dot) and the right junction (dot -drain lead) are related 
to each other, and the average occupation of the dot, by 

lL{t)-lR{t)^e'^. (24) 

In the steady state, the occupation of the dot is constant and the average currents through the two junctions are 
equal. Of course, the actual fluctuating time dependent currents are almost never equal. The external current arises 
as the external circuit adjusts the chemical potential of the local Fermi reservoir when electrons tunnel onto or off the 
dot. It is thus clear that the current through the left junction must depend only on the tunneling rates 7Li, i = 1, 2, 3. 
in the left barrier. This makes it easy to identify the average current through the left (or right) junction by inspection 
of the equation of motion for {c^c): all terms in the right hand side of Eq. (|24|l proportional to 7^4 correspond to the 
average current through the left junction, Ihit), while all terms on the right hand side proportional to ^Ri correspond 
to the negative of the current through the right junction, —Iii{t). 

IV. LOCAL SYSTEM DYNAMICS 

We can now compute the current through the quantum dot. The current reflects how the reservoirs of the source 
and drain respond to the dynamics of the vibrational and electronic degrees of freedom. Of course in an experiment 
the external current is typically all we have access to. However the master equation enables us to calculate the coupled 
dynamics of the vibrational and electronic degrees of freedom. Understanding this dynamics is crucial to explaining 
the observed features in the external current. As electrons tunnel on and off the dot, the oscillator experiences a force 
due to the electrostatic potential. While the force is conservative, the tunnel events are stochastic (in fact a Poisson 
process) and thus the excitation of the oscillator is stochastic. Furthermore the vibrational and electronic degrees of 
freedom become correlated through the dynamics. In this section we wish to investigate these features in some detail. 



From the master equation, the rate of change of this average electron number in the dot may be obtained: 

'-^ - ^^[^'^|] ^''^ 

= [7Li(l - X'KfiL (etc)) +7m(l - X'){hR ~ (etc)) 
-27LiA2(/iL(ata) - (a+actc)) 
+lL2\^{f2L{a^a) - (atactc) - (1 - /2L)(ctc)) 
+7L3A'(/3l(1 + ata) _ f^^(^c^c) ~ {a' ac^ c)) 
-27mA'(/iij(ata) - (atactc)) 
+7ii2A2(/2fl(ata) - (atactc) - (1 - /2/?)(ctc)) 

+7raA2(/3fl(l + flta) - /3fl(ctc) - (atactc))]cT • (26) 

While for the vibrational degrees of freedom, we see that 

'^(«ta)cT , r t "^^l ro7\ 

= A2[7i2(-/2L(ata) + (atactc) + (1 - /2L)(ctc)) 

+7L3(/3l(1 + at a) - /3L(ctc) - (atactc)) 

+7fl2(-/2fl(ata) + (atactc) + (1 - /2«)(ctc)) 

+7fl3(/3i?(l + ata) - /3_R(ctc) - (atactc)) - K(ata)]cT + Kfip . (28) 

Here the subscript CT indicates that the quantity to which it is attached is evaluated in the canonical transformed 
(CT) basis. The average occupational number of electron in the dot in the original basis is the same as in the CT 
basis: 

(etc) = tr[ctcp] — tr[ctcp] — (ctc)cT- 

While for the vibrational degrees of freedom, we have 

(ata) = tr[atap] 

= tr[atae-'^e-*"°''^''*pe''"'°''*''*e"] 

= tr[e"^°'^^''*(at + \n){a + \n)e-'^ •''''' '''■p\ 

= (ata)cT + A((ate-*"°'^'"* + ae-*"°'^''^*)n)cT + A2(n2). (29) 

If the initial displacement (x) is zero, the second time dependent term in the previous expression remains zero. We 
will assume this is the case in what follows. 

In general we do not get a closed set of equations for the mean phonon and electron number due to the presence 
in these equations of the higher order moment (atactc). This reflects the fact that the electron and vibrational 
degrees of freedom are correlated (and possibly entangled) through the dynamics. One might proceed by introducing 
a semiclassical factorization approximation by replacing (atactc) by the factorized average values, i.e., (atactc) w 
(ata) (etc), then the evolution equations (|26|l and (|28|) forms a closed set of equations. However there is a special case 
for which this is not necessary. If we let 7^1 — 7l2 = 7l3 — 1l and 7^1 = ^fa — 7_R3 — 1r which is the assumption 
of energy-independent tunnel couplings, the equations do form a closed set: 



and 



^^ - Ai(ctc)+i?i(ata)cT + Ci, 

Ai = -[il{1 - hL\^ + hL\^) +1r{1 - hnX^ + hR>?)\, 

Bl = X^{-2fiLlL + f2LlL + f3LlL-'2flRlR + f2R'lR + f3R.lR), 

Ci - (1 - X^hLfiL + IlIzlX" + (1 - X^)iRhR + IrHrX^ , 

d(ata)cT - , t X n / + V ^ 

^^ — = A2(ctc) + B2(ata)cT + C2 , 

A2 = A^(7L(l-/2L-/3L)+7fl.(l-/2i?-/3K)), 
B2 = X^{--lLf2L+lLfzL~lRf2R + lRhR)- H, 

C2 = X'^i'jLfsL + IrIzr) + i^np . 



(30) 



(31) 
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Consideration of Eq. (|31|l indicates that it is possible for the oscihator to achieve a steady state even when there is no 
exphcit thermo-mechanical damping (k = 0). This requires bias conditions such that f^^ = fan = 0. It is remarkable 
that the process of electrical transport between source and drain alone can induce damping of the mechanical motion. 
This result has been indicated by other authors^iiSiiiiiS. 

These equations were solved numerically and the results, for various values of A and bias voltage, are shown in 
Fig. 121 A feature of our approach is that we can directly calculate the dynamics of the local degrees of freedom, for 
example the mean electron occupation of the dot as well as the mean vibrational occupation number in the oscillator. 

From these equations we can reproduce behavior for the stationary current similar to that observed in the exper- 
iment. We concentrate here on the stationary current through the left junction (connected to the source). Similar 
results apply for the right junction. We assume that the electronic temperature is 1.5K, which is the temperature 
used in the experiment by Park et al^. 

Following the discussion below Eq. (|24|l . we see from Eq. H3U|) that the average steady state current through the 
left junction is given by 

1st = e7L[(-l + /2lA2 - /3LA2)(ctc},t + (-2/ilA2 + f^^X'' + f3LX^){a^a)cT,st 

+ (l-A2)/ii + /3LA2]. (32) 

which is of course equal to the average steady state current through the right junction. The steady state current, 1st 
can then be found by finding the steady state solution for each of the phonon number and electron number, 

+ , B1C2 - B2C1 

^''^^' ^ A,B2 - A2B, ' (^3) 

, t \ ~A2 f B1C2 - B2Ci\ C2 , . 

^" °^^^'^' ^ -W Ub2^A2bJ - B-2 • (^^) 

In Fig. 121 we assume that the bias voltage is applied symmetrically, i.e, fiL — —fJ-R — eVbias/2. In this case, 
all the Fermi factors fm, f2R, and fsR effectively equal to zero in the positive bias regime, the regime of Fig. |21 
From Eqs. (|24|l and (|30|l . we see that Eq. H32I) also equals to the steady state current through the right junction as: 
^st = e7fl(c^c)st- In the case 7/j = 7^,, the steady state (c^c)st shown in Fig. I^Ja) and (d) at long times should thus 
equal respectively to Ist/{&lL) shown in Fig.[2Ic) and (f) at long times. This is indeed the case, although the transient 
behaviors in these plots differ considerably. We note that the plot shown in Fig. |2Ic) and (f ) is the current through 
the left junction, normalized by {e^l)- The values of /il, /2L, and /3L depend on the strength of the applied bias 
voltage and are important in understanding the stepwise behavior of the stationary current as a function of the bias 
voltage. We will now concentrate exclusively on the positive bias regime. 

When the bias voltage is small, the current is zero. As the bias voltage increases the first Fermi factor in the left lead 
to become non zero is /2L, with the other Fermi factors very small or zero. In the case that /2L = 1, /il = 0, /ai = 0, 
the steady state current is 

4^) = e7i[A2(ata)cT.st - (1 - \^){c^c)st] , (35) 

For low temperatures this is very small. Only if the phonon temperature is large, so that the stationary mean phonon 
number is significant, does this first current step become apparent (see Fig. 0J|. As the bias voltage is increased both 
/2L and fiL become non zero. In the case where they are both unity, the steady state current is 

/i') = e7i[(l - \^){cc^U - A2(ata)cT,st] • (36) 

The first term here is the same form as the bare tunneling case except that the effective tunneling rate is reduced by 
(1 — A^). This is not too surprising. If the island is moving, on average it reduces the effective tunneling rate across 
the two barriers. Thus the first current step will be reduced below the value of the bare (no phonon) rate. At the 
region where bias voltage is large, all the Fermi factors are unity and 

/if) = e7L(cct),, , (37) 

which is the expected result for the bare tunneling case. 

To explicitly evaluate the stationary current we need to solve for the stationary mean electronic and phonon 
occupation numbers. We have done this numerically and the results are shown in the figures below. However the 



large bias case can be easily solved; 



IL 



{c^cU = ^^— , (38) 

{a^aU = U + ^!iz7£±M) {^11— \ + 1I^^L±^ , (39) 

1st = e ^ . (40) 

IL +1R 

This is the result for tunneling through a single quasi bound state between two barriers^^. 

The steady state current for larger values of A shows two steps. As one can see from Fig.EJa), the current vanishes 
until the first Coulomb Blockade energy is overcome. The first step in the stationary current is thus due to bare 
tunneling though the dot. The second step represents single phonon mediated tunneling through the dot. These 
results are consistent with the semiclassical theory of Boese and Schoeller— given that our expansion to second order 
in A can only account for single phonon events. The height of the step depends on A which is the ratio of the coupling 
strength between the electron and the vibrational level x, and the oscillator energy htOo- 

Looking at Fig.|21 the average electron number approaches a steady state (e.g., a steady state value of 0.5 at large 
bias since we have set the value of 7l to be equal to 7^) while the average phonon number, without external damping, 
behaves differently in various regions (Fig. [Sfb)). The average phonon number slowly reaches a steady state value 
within the first step, while at the second step, the mean phonon number grows linearly with time (Fig. |2Ib)). The 
steady state at the first step is the previously noted effect of transport induced damping. These results are as expected 
since from Eq. (|2f|l . the term consisting ^l2^'^ f2L'D[ac']p corresponds to a jump of one electron onto the island with 
the simultaneous annihilation of a phonon, while 7L3A^/3iX'[a^^c^]p corresponds to the jump of an electron onto the 
island along with the simultaneous creation of a phonon. The dynamics can be understood by relating the behavior 
of these terms to the rate of average phonon number change in Eq. (|31l) . At the first step, when both /i^ and /2L 
are both one, while f^L is zero, the coefficient B2 has negative value, and therefore the mean phonon number could 
reach a steady state under this transport induced damping. In this regime, we find that 

(ctc).t = (1-A2)/2, (41) 

{a^aU « 1/2, (42) 

where we have set jl — 1r and used Eq. l(^ in obtaining Eq. (jl^ . The corresponding effective temperature can be 
found using 

^'"^ " fcsln[l + (l/(ata),0] ' ^^^^ 

When all the Fermi factors for the left lead are unity, the rate of growth for the mean phonon number now depends on 
a constant C2 in Eq. (|31|1 and therefore the mean phonon number will grow linearly with time. However, the current 
will be still steady (Fig. I3f)). This indicates that the steady state current and mean electron number in the dot 
[given by Eqs. H40|l and H38|) respectively] do not depend on the mean phonon number. This is supported by the fact 
that the coefficient Bi in Eq. (|30|) vanishes in this regime. When damping is included, the phonon number reaches a 
steady value of 0.35 [see Fig.|5fe) and also Eq. (|^ ]. 

In Fig. |31 we plot the steady state current versus bias voltage for different values of A. At the region of the first 
phonon excitation level (when bias voltage is between 30 meV and 40 meV), the steady state current drops by a factor 
proportional to quadratic order of A. Compared to the current at large bias voltage, the size of the drop is 

.r _ 7LA^(27LA^ + ^(2np + l)) 

^^^* - 47lA2-2^(A2-2) • ^^^) 

We thus see that the effect of the oscillatory motion of the island is two fold. Firstly the vibrational motion leads 
to an effective reduction in the tunneling rate by an amount A^ to lowest order in A. Secondly there is a second step 
at higher bias voltage due to phonon mediated tunneling. This is determined by the dependence of the Fermi factors 
on the vibrational quantum of energy. One might have expected another step at a smaller bias voltage. However this 
step is very small unless there is a significant thermally excited mean phonon number present in the steady state. If 
we increase the phonon temperature such that the energy is larger than the energy quantum of the oscillator, hujo, (for 
this example we choose T = 2huJo/kB ~ 116i^ ) this step can be seen (Fig^. Thus we see three steps as expected, 
corresponding to the three bias voltages when the three Fermi factors switch from zero. 

In order to explore the steady state correlation between phonon number and electron number on the dot, we can 
find the steady state directly by solving the master equation from Eq. H21|l . In Fig. Owe plot the correlation function 
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FIG. 2: Average number of electron, phonon and current through the dot against bias vohage with hujo ~ 5 meV, hui — rj 
= 15 meV, ksT = 0.13 meV, and h'yL = h'yR = 2 /leV for A = 0.3. Figures (a),(b),(c) without damping and (d),(e),(f) with 
damping k = 0.37l. We assume fj,L — —jj-r — eVbias/2. 

(a'ac'c) — {a'a){c'c) as a function of A an(i the bias voltage. The correlation is seen to be small except when a 
transition occurs between two coniiuctance states. This is not surprising, as at this point one expects fluctuations in 
the charge on the (iot, and consequently the fluctuations of phonon number, to be large. This interpretation implies 
that damping of the oscillator should suppress the correlation, as the response of the oscillator to fluctuations in the 
dot occupation are suppressed. This is seen in Fig. O 



CONCLUSIONS 



We have given a quantum description of a QEMS comprising a single quantum dot harmonically bound between two 
electrodes based on a quantum master equation for the density operator of the electronic and vibrational degrees of 
freedom. The description thus incorporates the dynamics of both diagonal (population) and off-diagonal (coherence) 
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FIG. 3: Steady state current for different values of A with damping k = 0.37l, the electronic and phonon temperatures are 
both 1.5K. 
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FIG. 4: Steady state current for different values of A and damping k = O.371,, the electronic temperature is 1.5K, and the 
phonon temperature is 116K which is chosen to be 2hjjJo/kB in order to make manifest the step at a smaller bias voltage. 



terms. We found a special set of parameters for which the equations of motion for the mean phonon number and the 
electron number form a closed set. From this we have been able to reproduce the central quahtative features of the 
current vs. bias voltage curve obtained experimentally by Park et al.^ and also of the semiclassical phenomenological 
theory by Boese and Schoeller^^. We also calculate the correlation function between phonon and electron number 
in the steady state and find that it is only significant at the steps of the the steady state conductance. The results 
reported in this paper do not probe the full power of the master equation approach as the model does not couple the 
diagonal and off-diagonal elements of the density matrix. This can arise when the vibrational motion of the dot is 
subject to a conservative driving force, in addition to the stochastic driving that arises when electrons tunnel on and 
off the dot in a static electric field. The full quantum treatment will enable us to include coherent effects which are 
likely to arise when a spin doped quantum dot is used in a static or RF external magnetic field. 
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FIG. 5: Difference in the correlation function for different values of A with damping k = 0.37l- 
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FIG. 6: Difference in the correlation function for different values of k with damping A 



0.37L. The plot starts at a value 



Acknowledgments 

G.J.M and D.W.U acknowledge the support of the Australian Research Council Federation Fellowship Grant. H.S.G 
would like to acknowledge support from Hewlett-Packard. We gratefully acknowledge discussions with Jason Twamley, 
supported by EC 1ST FET project IST-2001-37150 QIPDF-ROSES. 



^ M. Roukes, Physics World 14, 25, February (2001). 
2 R.G. Knobel and A.N. Cleland, Nature 424, 291 (2003). 

^ M. D. LaHaye, O. Buu, B. Camarota, and K. C. Schwab, Science 304, 74 (2004). 

* Z. Zhang, M. L. Roukes, and P. C. Hammel, Journ. App. Phys. 80, 6931 (1996); H. J. Mamin, R. Budakian, B. W. Chui 
and D. Rugar, Phys. Rev. Lett. 91, 207604 (2003). 



13 

T. A. Bmn and H.-S. Goan, Phys. Rev. A 68, 032301 (2003); G. P. Berman, F. Borgonovi, H.-S. Goan, S. A. Gurvitz and 

V. I. Tsifrinovich, Phys. Rev. B 67, 094425 (2003). 

A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. Rev. Lett. 88, 148301 (2002). 

X. M. H. Huang, C. A. Zorman, M. Mehregany, and M. Roukes, Nature 421, 496 (2003). 

H. Park, J. Park, A. K. Lim, E. H. Anderson, A. P. AUvisator, and P. L. McEuen, Nature 407, 57, (2000). 

D. Mozyrsky and I. Martin, Phys. Rev. Lett. 89, 018301 (2002). 

A. Yu. Smirnov, L. G. Mourokh, and N. J. M. Horing, Phys. Rev. B 67, 115312(2003) 

D. Mozyrsky, L Martin and M. B. Hastings, Phys. Rev. Lett. 92, 018303 (2004). 

A. D. Armour, M. P. Blencowe, and Y. Zhang, Phys. Rev. B 69, 125313 (2004). 

D. Boese and H. Schoeller, Europhys. Lett. 54, 668 (2001). 

H. B. Sun and G. J. Milburn, Phys. Rev. B. 59, 10748 (1999); G. J. Milburn, Aust. J. Phys. 53, 477 (2000). 

V. Aji, J. E. Moore and C. M. Varma, Electronic-vibrational coupling in smgle-molecule devices, cond-mat/0302222 (2003). 

C. W. Gardiner and P. Zoller, Quantum Noise, 2nd edition, (Springer- Verlag, Berlin, 2000) 

R. I. Shekhter, Yu. Galperin, L. Y. Gorelik, A. Isacsson and M. Jonson, J. Phys.: Condens. Matter 15 R441 (2003). 

G. Mahan, Many-Particle Physics, (Plenum Press, 1990). 



